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Abstract. A new calculation of -RdAu has been performed using the 2003 d-|-Au data and the 
higher-statistics 2005 p-l-p data. These nuclear modification factors are compared to calculations 

fljl ' using nuclear-modified PDFs and a J/i/' breakup cross section is extracted. These values are 

then used to project the cold nuclear matter efi^ects in Au+Au collisions. Additionally, a more 

fj . data-driven projection is performed. 

^ ■ 1. Introduction 

^D ! Cold nuclear matter (CNM) effects are those due to bemg within the nuclear environment, as 

^^ ' opposed to a higher-density medium or vacuum. These effects are quite interesting on their own, 

,—1. . but are even more important at the Relativistic Heavy Ion Collider (RHIC) as we attempt to 

observe the Quark-Gluon Plasma (QGP), a state with much higher energy density than normal 
(^ ■ nuclei. There are several proposed signatures that could indicate the transition to a QGP, and 

OO ! it is essential to understand them in the context of what normally occurs in nuclei. 

^^ ' One such proposed signature is the suppression of J/V' production in the deconfined medium 

of the QGP, due to color-charge screening of the cc interaction [T] . However, J /ip production 
k>( \ may also be modified due to CNM effects such as nuclear shadowing, gluon saturation, the EMC 

\^ • effect, or absorption/breakup within the nuclear remnants. Therefore it is very important to 

^ quantify these contributions. 

2. The PHENIX Experiment 

In order to study CNM effects at RHIC, deuteron and gold ions are collided at ^/snn = 200 
GeV. This was first done in Run-3 and more recently in Run-8, the latter of which is currently 
being analyzed. Since the original analysis of the Run-3 data [2], a new p-|-p dataset was 
recorded in Run-5 with more than an order-of-magnitude increase in J/ip statistics over the 
p-|-p data from Run-3 [3]. In addition, there have been two years' worth of improvements in 
the reconstruction software, the J/'0 signal extraction, and our understanding of the detector. 
Consequently, a new analysis of the Run-3 d-|-Au data was performed using the same methods 
as the new p-|-p reference data so that an apples-to-apples comparison could be done [1]. It 
should also be noted that the Run-4 Au-|-Au [5J and Run-5 Cu-|-Cu 6J Raa results use this 
same p-|-p reference data. 



The PHENIX rapidity coverage consists mainly of two regions: the Drift Chamber, Pad 
Chamber, EM Calorimeter and RICH measure J/ij: -^ e^e~ at mid-rapidity {\y\ < 0.35), and 
the Muon Tracker and Muon Identifier measure J/'0 -^ A^^/^~ ^^ forward and backward rapidity 
(1.2 < \y\ <2.2). 



3. Nuclear Modification Factor 

The nuclear modification factor .Raa (Equation [1]) quantifies the suppression or enhancement 
of particle production in collisions of heavier nuclei with respect to p+p collisions, scaled by 
the appropriate number of binary collisions ((A'coii)) in the heavier species, as calculated by a 
Glauber model [7]. Using the Run-3 d+Au and Run-5 p+p data we can construct -RdAu as 
functions of transverse momentum, collision centrality, and rapidity. We will focus on the latter 
here; for the broader analysis, see ^. 
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-RdAu as a function of rapidity is plotted as the black points in Figure [TJ As can be seen in 
the Figure, J/'0 production in d+Au collisions is consistent (within the large error bars) with 
binary-collision-scaled p+p collisions at backward and mid-rapidity, but is suppressed at forward 
rapidity (in the gold-going direction). 




1 


- 


T 


m 






• 


:--| 






i 


-I 











11% Global Scale Uncertainty 



< 

■D 



0.5 




NDSG Model 



breakup 



Best Fit a,,,,,„p = 2.8 :;>b 



breakup 



: 0,1 ,2,3,5 mb (top to bottom) 



Best Fit a,,,,,,p = 2.2 :;|mb 



-2 



-1 

Rapidity 



-2 



-1 
Rapidity 



Figure 1. (color online) RdAu data compared to various theoretical curves for different (Threakup 
values. Also, shown as a band are the range of (Jbreakup found to be consistent with the data 
within one standard deviation. The left panel is a comparison for EKS shadowing [8], while the 
right panel is for NDSG shadowing [9]. 



It can be useful to compare the measured i?dAu to a simple model of CNM effects. By 
comparing the modification due to available nuclear-modified PDFs with the data, we can 
extract a breakup cross section {(Jhreakup) for the J/ip passing through the nuclear medium. 
This is found to be Chreakup = '^■^-ia ™b and abreakup = 2.2;1^^;5 mb using the EKS [8] and 
NDSG [9] nuclear-modified PDFs, respectively. Note that the quoted uncertainties account 
for all experimental statistical and systematic errors, including the global scale uncertainty. 



However, theoretical model uncertainies are not included. For details of how the statistical and 
systematic errors of the data are accurately accounted for, see [1] and [TO] . 

These values overlap within one standard deviation with the published value of cr^j^J = 
4.5 lb 0.5 mb from CERN-SPS [IT] (within the large uncertainties). It should be noted, however, 
that the SPS value does not include any nuclear shadowing or anti-shadowing. It has been 
suggested that the SPS data lies in the anti-shadowing region, causing an enhancement in J/ip 
production, which would then have to be balanced by an even larger 0"^^^, to match the data [|12j. 

4. Projections to Au+Au 

To quantify CNM effects in Au+Au collisions, we use the nuclear-modified PDFs (in conjunction 
with a model of their impact-parameter dependence [l3l E]) and the calculated (Jbreakup to 
project -Raa for Au+Au strictly due to CNM effects. This is shown versus the number of 
participant nucleons in the Au+Au collision in Figure O where the two bands represent the 
one-standard-deviation uncertainty bands using the two PDFs, and the blue points represent 
the PHENIX J/ip Raa measurements from Run-4 [5]. 

The -Raa measured in Run-4 shows statistically-significant suppression at forward rapidity 
in comparison to the projected suppression due to CNM effects. This is not true, however, in 
the mid-rapidity case, except perhaps in the most central collisions. There are several things to 
note here: first, that the large error bands on the CNM projections limit our ability to make 
any definitive statements, and this will have to be improved in the future; second, that the 
above calculations are explicitly model-dependent, and there is no easy way to include this in 
the uncertainty bands; third, the error bands shown in Figure [5] are correlated between the two 
rapidity regions, as they are due to the same calculation of criyreakup- Calculations in separate 
rapidity bins were performed for [3], but are not included here due to the large uncertainties. 
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Figure 2. (color online) Raa for Au+Au [5] collisions compared to a band of theoretical curves 
for the CFhreakup values found to be consistent with the d + Au data as shown in Figure [H The 
left figure includes both EKS shadowing [8j and NDSG shadowing [9j at mid-rapidity. The right 
figure is the same at forward rapidity. 



A less model-dependent approach to the problem is to parametrize the CNM data directly, 
as proposed in |15j . Instead of invoking nuclear-modified PDFs and breakup cross sections, we 
start from the assumption that i?dAu should approach unity for the most peripheral collisions, 
and then fit the measured i?dAu fo a linear function satisfying this condition (convolving this 



function with the d+Au centrality distributions for each data point). Other fit functions may 
be tried, but the results do not vary much due to the imprecision of the current data. 

Using the results of the fit we then construct .Raa as a function of impact parameter by 
integrating over the radial distribution of binary collisions for each gold nuclei at a given impact 
parameter value. As can be seen in Figure [3l the resulting one-standard-deviation bands are 
qualitatively similar in shape to those shown in Figure EJ but generally have a larger la region 
and slightly more suppressed most-likely values. It must be noted that, contrary to the previous 
case, the uncertainty bands here are entirely uncorrelated between rapidities. 
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Figure 3. Predictions of the data-driven method constrained by the -RdAu a-s a function of 
collision centrality for the Au-|-Au Raa for mid-rapidity (left) and at forward rapidity (right). 



5. Summary 

We have presented RdAu as a function of rapidity as a quantification of CNM effects in heavy 
ion collisions. We have also used several techniques to project the CNM effects to -Raa for 
Au-|-Au collisions. 

MGW acknowledges funding from the Division of Nuclear Physics of the U.S. Department of 
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